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Abstract 

A new analytical criterion that captures the onset of separation of flow past elliptic cylinders 
is derived by considering the variation of the wall normal velocity in Reynolds number parameter 
space. It is shown that this criterion can be used to calculate the separation Reynolds number 
(Rcs) for the classical problem of flow past a circular cylinder, a contentious and unresolved issue 
till date. The two dimensional Navier-Stokes equations are solved computationally and an exact 
value of Rcs is obtained by applying the aforementioned criterion. 
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A great deal of research in the past century has been focused on bluff body wakes and 
considerable progress has been made towards understanding them. A body is considered 
"bluff' if the spatial extent of the body along the flow direction is of comparable or lesser 
order to that normal to it. The bluff bodies that are of concern here are a family of symmetric 
elliptical cylinders, with their minor axes aligned with the flow direction. The flow past such 
a bluff body is steady for very low Reynolds numbers {Re = Ua/u where a is the length of 
the body normal to the flow and u is the coefficient of viscosity), which in case of a circular 
cylinder happens for Re < 47. In two dimensions, as Re is increased, a few well defined 
features of the flow are observed. Up to a particular value of Re (say Res), the streamlines, 
while being asymmetric about the axis normal to that of the flow, are attached to the 
body 1]. But for Re > Reg, the flow separates and two well defined separation "bubbles" 
or eddies are are regions of closed streamlines where the 

flow direction just next to the cylinder is in the opposite direction to the mean flow {i.e. 
a region of backflow). The characteristic features of these bubbles have been documented 
comprehensively both in computational and experimental studies, mostly for the case of flow 



past a circular cylinder. Investigations of flow past elliptic cylinders are few J2] and none of 
them seem to focus on aspects of flow separation. Computational results ^, IJ indicate that 
'or the case of a circular cylinder, 4 < Reg < 7, but even recent investigations, for example 
511, have been unable to capture the bubble for Re < 6 inspite of having a much higher 
„„.e*a. .solution in tKe wa.e .e,on than S Q. Expe..eMa> .udies fl Q su.e. 
from similar ambiguities and an effort to summarize and compile various such experiments 
was undertaken by [9] who concluded that Re.^ ~ 5. The criterion derived here, as we shall 
see later, can be used to obtain the exact value of Res and resolve this ambiguity. 

The complexity of the problem means that only a few comprehensive theories have been 
put forth in the past and their scope and success have been limited. One theory of note 



by [Id, ini ll3 uses the triple-deck model to show that the length of the bubbles (in the 



his theoretical result has 
However, this theory 



direction of the flow) increases linearly with Reynolds number, 
been confirmed by many experiments and numerical simulations 
is unable to predict the Reynolds number at which the bubbles appear. Res- Further, it does 
not attempt to explicitly understand and characterize the flow when bubbles form. Neither 
does it attempt to generalize the same for other bluff bodies like symmetric elliptic cylinders 
in which the bubble formation is of a similar nature. 
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FIG. 1: Schematic of the flow geometry 

The aim of this Rapid Communication is an attempt, a possible first step, in this direction. 
Here we derive a simphstic analytical criterion that captures the Reynolds number at which 
the bubbles start forming, i.e. Rcg. This criterion is general and valid for a whole class of 
elliptic cyhnders. It seems similar to the Prandtl shear stress criterion which states that, 
at the point where the streamline separates the wall shear stress, r^aii = 0. The distinction 
between the Prandtl criterion and the present one though is quite fundamental and should 
pose no source of confusion. The criterion derived here is valid only at Rcg, i.e. it captures 
the Reynolds number at which the bubble starts forming. Whereas once the bubble forms, 
the upper and lower streamlines ending at the wall (the "separatrices" since they separate the 
"outer" flow region from the flow "inside" the bubble) always satisfy the Prandtl shear stress 
criterion at the point that they separate, irrespective of the Reynolds number (of course, 
until the onset of unsteadiness). The Prandtl criterion, it must be noted, has relevance only 
after the bubble is formed and the point that the separatrices separate from the cylinder 
surface has Tu;aii — 0, but this has absolutely no bearing on Rss. 

We study uniform flow past a 2D elliptic cyhnder such that the flow direction is along 
the minor axis, as shown in Fig. 1. A body-fltting orthogonal coordinate system given by 
C — 1] (specifically, an elliptic cylindrical coordinate system) coincides with the x — y axis 
at the point P, the intersection of the symmetry plane with the cylinder. One expects that, 
as the Re is increased from zero through Res, it is at P that the bubbles start forming. 
Examining the flow in the vicinity of P in Re parameter space would give us the required 
criterion. This is done by considering the flow around P for two Reynolds numbers Rei and 
Re2 such that 

< Rci < Rcs < Re2 < Rcuns (1) 
Here Reuns is the Reynolds number at which the steady wake of the cylinder becomes 
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unstable and vortex shedding sets in, e.g. Rcuns — 47 for a circular cylinder. Now for 
Re — Rei, the bubble has not started forming and there is no region of backflow; whereas 
for Re = Re2, the bubble has formed and is of finite size and has a region of backflow. Then 
one can find a distance which is less than the bubble length, so that. 



u{C,p + K,VP,Rei) > 



(2) 
(3) 



Here P = {Cp,i]p), P' = {Cp + 6(,rip) and the flow velocity in ( — rj coordinates is written 
as M = -u^ and v = Ur, along the respective coordinates. 

Note that the no slip condition at the wall implies that Vl\p—y\p—Q. Since v — Q along 
the cyhnder, |^|p = (so are derivatives to higher orders along that direction). Using the 
continuity equation at P, ^§^ + ^f| — (here /i^ and /i^ are the scale factors for the 



coordinate transformation (x,y) 



(C, ry) ) and the fact that |^|p = 0, we get 



du 



= 



(4) 



Now we write the u velocity at point P' as a Taylor expansion about the point P, keeping 
terms till second order. Since 5( is along the ( direction, we get 



U\p' 



. du 



Prom no shp and Eq. (4) we get 



U\pt 



, 2! 



(5) 



(6) 



Clearly, the relations (2) and (3) imply that u\pi{Re)—0 at some Rei < Re < Re2. By 
letting the hmits Rei — > ReJ and Re2 Ref alongwith S( 0, we see that this Reynolds 
number is precisely Reg, i.e. u\p'{Res)—0. And from (6) we finally get. 



dC 



(7) 



which is precisely the criterion required for separation. The above derivation is almost trivial 
and the result is fairly intuitive. Purther, it is seen that the condition (7) is necessary and 
sufficient. This is because u near P is always positive for Re < Res and so is the second 
derivative. Similarly for Reuns > Re > Reg it is always negative, again implying a negative 



value of the second derivative. An equivalent pressure condition can also be derived by 
using (7) and considering the Navier-Stokes (NS from here on) equation at P. Setting all 
derivatives of velocity components along the wall to be zero and using Eq. (4), we get, 

1 d'^u dp 

(Cp, Vp, Re) = ^e— (Cp, r/p, Re) (8) 



for Re = Rcg, we have from (7), 

|^(Cp,r/p,i?e,)=0 (9) 

Of the two equivalent conditions, (7) and (9), it is not clear which would be more useful in 
computational efforts in determining Reg- This can be resolved by examining the behavior 
of ^{Cp,r]p, Re) and ^{Cp,r]p, Re) in Re space around Res- 
It is instructive to note that relations similar to (7) and (9) would also be valid in 
three dimensions for flow past a family of symmetric ellipsoids. This family, to which the 
sphere also belongs, has rotation symmetry about an axis parallel to the flow direction. 
The derivation of the separation criteria for this family mirrors the one above. While flow 
past an ellipsoid has not been investigated in any detail, for a sphere, experiments imply 
that Res=25 jl^ . 15^ . However, unlike the case of a flow past a cylinder, there is excellent 
agreement between various experimental and computational results jl^ about this value. 
In the light of this fact, the criterion (7) has little application here. Our present study is 
therefore restricted to the study of two dimensional flows. 

We consider flow past a circular cylinder, which is a special case of the family of ellipses 
and relatively easy to compute. In this case, {(^v) = ^) but we continue to use {(^v) 
for the sake of uniformity of notation. We solve the 2D NS equations computationally in 
the streamfunction-vorticity formulation. This approach has an advantage over primitive 
variable formulations since an explicit operator splitting is not required. The unsteady NS 
equations are solved in the range 2.5 < Re < 10 by impulsively starting the flow, unlike 
the authors 4] who solve the steady NS equation. Impulsively started flow past a cylinder 
has been studied extensively using a variety of numerical formulations - flnite difference, 
flnite volume and vortex methods [la . The present work uses a flnite difference formulation 

n 

for its simplicity as done by p||. It is formulated using an explicit time stepping scheme, 
the time stepping being varied from a fourth Runge Kutta scheme to an Euler scheme 
with identical results. The discretization of the non-linear terms was done using a 3rd 
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distance from P 

FIG. 2: Variation of the u velocity along the symmetry line starting from P at different Reynolds 
Numbers. 

order upwind scheme while all other derivatives were effected using central differencing. 
The streamfunction-vorticity Poisson equation was solved using a stabilized Biconjugate 
Gradient Method [l^. The grid used was a body fitting grid, with grid clustering in the 
radial direction. Computations were performed on three different grids with the highest 
resolution being 300 x 150 and the largest ratio of the outer boundary to the cylinder 
diameter being 40. Each computation was performed till steady state was reached. About a 
hundred hours of computation time were required for all the cases on the most refined grid, 
the results of which are presented below. 

Once the flow is started impulsively from rest, a pair of separation bubbles can be seen 
for Re > 6, below which it was hard to resolve the bubbles. The computations were 
continued up to a time t=10 (where U = = 1, U and D being the free stream velocity 
and cylinder diameter respectively) when steady state was approximately reached. The 
difficulty in resolving the bubbles close to separation is not a new problem and is common 
to computational investigations of this nature In order to do this, we first look at the u 
velocity variation along the symmetry plane close to the cylinder. This has been done and 
the results are plotted in Fig. 2 for a range of Reynolds numbers as indicated. Clearly u < 
near P indicates the presence of a bubble. As seen from Fig. 1, 4 < Rcs < 6, which is in 
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FIG. 3: Variation of |^(Cp, ??p, i?e) with Reynolds number at the point P. 

good agreement with existing computational results. 

In Fig. 3, ^(Cp, VP: -^c), which is obtained from the u velocity data by central differenc- 
ing, is plotted against Reynolds number. As seen in Fig. 3 we obtain the rather surprising 
result that ^{(p,r]p, Re) varies linearly with Re in the neighborhood of Reg and so, by 
virtue of (8), ^{(p,r]p, Re) does not. Therefore, numerical simulations attempting to com- 
pute Res can do so by calculating ^{(p,r}p, Re) for a few Reynolds numbers around Reg 
and then extrapolating the curve linearly to zero to obtain Reg- §^(Cp) Vp^ R^) ^^e other 
hand is highly non-linear around Reg and cannot be used in similar manner, even though (9) 
itself remains valid. So from Fig. 3 and applying (7) one gets it!es=4.8. The aforementioned 
approach underlines the advantage of using (7) in order to get Reg as opposed to attempts 
in resolving the bubble by solving the NS equation for a large number of values of Re. One 
must realize that, though (7) is exact, the linearity of |^(Cp, Vp-, R^) around Reg is observed 
only through computations and that too for the special case of a circular cylinder. But from 
the seemingly general nature of the problem of separation for the family of ellipses, as seen 
in deriving (7), one might expect that the linearity property holds well across the family. 
But an ad hoc argument of generality can hardly suffice and therefore the author is currently 
investigating this aspect of the problem, both analytically and numerically. 

Also, a note of caution regarding the value of Reg obtained here, which it must be re- 
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marked is still contingent on the computational approach used. And the present computation 
while reasonable, can be improved with respect to the grid refinement and also the outer 
domain size, which might have some bearing on the value of Rss- In fact there was some 
dependence of the exact value of Rcg on the outer domain size (about 7% between 30 and 40 
times the cylinder diameter for the most refined grid). But the linearity property remained 
unchanged and the variation in the slope of the ^(Cp, ^p, -Re) curve was also negligible. 
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